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ABSTRACT 

Preliminary  results  of  an  experiment  to  determine  the 
accuracy  and  the  effects  of  varying  amounts  of  environmental 
data  at  the  NAVTORPSTA,  Keyport ,  Washington,  Dabob  Bay 
Range  are  presented.   A  75-KHz  source  was  tracked  at  six 
depths  at  each  of  six  horizontal  ranges.   The  data  were 
analyzed  by  the  NAVTORPSTA  NUTRACK  III  isovelocity-layer 
computer  program  and  by  an  isogradient-layer  program 
developed  here,  STUTRACK  I.   The  results  show  that  a  low 
signal-to-noise  ratio  is  not  a  problem  even  at  long  distances 
and  that  an  isogradient  analysis  yields  superior  results 
with  less  environmental  data  than  does  the  isovelocity 
method. 
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I.   INTRODUCTION 

A.   DABOB  BAY  FACILITY 

A  facility  for  the  three-dimensional  tracking  of 
surface  ships,  submarines,  and  torpedoes  is  located  at 
the  Dabob  Bay  facility  of  the  Naval  Torpedo  Station 
(NAVTORPSTA) ,  Keyport,  Washington.   Chosen  because  of  its 
favorable  oceanographic  features  and  its  proximity  to  the 
NAVTORPSTA,  Keyport,  Dabob  Bay  has  been  used  for  torpedo 
testing  since  19^9.   It  consists  of  a  250-KHz  "high 
frequency"  installation  no  longer  operating,  and  a  75-KHz 
"low  frequency"  primary  tracking  system  enveloping  nearly 
the  entire  bay  providing  for  underwater  tracking  in  a 
volume  approximately  30,000  feet  by  4,500  feet  by  600  feet 
deep. 

Dabob  Bay,  a  branch  of  Hood  Canal,  is  in  a  deep  (600  feet) 
depression  adjacent  to  the  Olympic  Mountains.   The  Quilcene 
River  flows  into  the  northernmost  part  and  the  Dosewallips 
River  into  the  southernmost.   The  sides  of  the  bay  are 
precipitous  and  predominately  rocky,  while  the  bottom  is 
mud.   Annual  precipitation  averages  about  51  inches  and 
while  average  snowfall  is  about  15  inches,  any  one  snowfall 
seldom  remains  longer  than  a  few  days.   Winds  generally  blow 
along  the  length  of  Dabob  Bay;  i.e.,  North  and  South. 
Southerly  winds  are  usually  accompanied  by  clear  skies  and 
occur  most  often  in  summer.   Normal  wind  velocities  are 
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5  to  15  miles  per  hour.   Tidal  levels  in  Dabob  Bay  range 
from  about  -5  to  +15  feet. 

Water  temperatures  change  markedly  from  season  to  season 
and  even  from  day  to  day.   The  changes  are  most  noticeable 
near  the  surface  but  may  be  observed  even  at  greater  depths. 
Salinity  varies  considerably  with  depth  as  fresher  surface 
water  from  terrestrial  runoff  overlays  more  saline  water. 
The  variation  is  also  quite  seasonal,  reaching  a  maximum  in 
the  spring.   The  main  contributor  to  sound  velocity  changes 
is  water  temperature,  and  the  resulting  sound  velocity 
profile  is  representative  of  the  temperature  gradient. 

B.   MOTIVATION 

NAVTORPSTA,  Keyport ,  has  developed  a  computer  program, 
NUTRACK  III,  which  divides  the  bay  into  laterally  homogeneous, 
isovelocity  layers  and  traces  a  sound  ray  from  a  three- 
dimensional  hydrophone  array  mounted  just  above  the  bottom, 
back  to  its  source. 

An  experiment  was  undertaken  to  provide  comparison  of 
the  three-dimensional  position  of  a  source  with  respect  to 
a  receiver  as  determined  by  various  ray-tracing  techniques. 
From  this  comparison  it  was  hoped  that  an  optimum  ray-trace 
procedure  would  be  indicated  and  that  a  better  understanding 
of  the  significance  of  environmental  variations  would  result. 
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II.   PROCEDURE 

A.   DETAILED  RANGE  DESCRIPTION 

The  basic  underwater  tracking  system  on  the  75-KHz 
range  is  comprised  of  three  components:   hydrophone  array, 
computer  system,  and  vehicular  instrumentation.   The  ship- 
board transducer  emits  a  75-KHz  acoustic  pulse  in  synchronism 
with  a  master  clock  at  the  computer  site.   This  pulse  is 
detected  by  each  of  the  four  hydrophones  (RX>R  ,R  and  R  ) 
in  the  array  and  mixed  with  a  local  oscillator  having  a 
frequency  unique  to  that  hydrophone.   The  outputs  of  these 
oscillators  are  preamplified  and  fed  to  a  multiplexer,  which 
sums  the  signals  and  transmits  the  composite  signal  via  an 
underwater  cable  to  the  computer  site,  where  it  is  processed 
to  determine  the  in-water  transit  time  to  each  hydrophone 
(Tx,Ty,Tz,  andTc). 

1.   Hydrophone  Array  Configuration 

The  hydrophones  in  the  NAVTORPSTA  Range  are  arranged 
in  a  short-baseline  system  in  groups  of  four,  each  one 
located  on  each  of  four  adjacent  corners  in  an  imaginary 
30-foot  cube,  thereby  defining  the  orthogonal  coordinate 
system  in  which  the  measurements  were  made.   These  hydrophone 
arrays  used  buoyant  spheres  which  exert  an  upward  force  in 
excess  of  two  thousand  pounds  to  keep  them  upright.   The 
array  is  not  free  to  rotate  about  the  Z  axis,  but  may  tilt 
as  much  as  two  degrees  from  the  horizontal  in  strong  currents. 
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This  tilt  is  measured  by  servo-pendulum  transducers  whose 
specified  resolution  is  O.OO833  degrees,  accuracy  0.025 
degrees,  and  linearity  0.05  percent  over  the  full  scale. 

The  Range  consists  of  six  arrays  (Pig.  1)  aligned 
along  a  line  ^50  feet  east  and  parallel  to  the  Range  axis, 
spaced  at  6000-foot  intervals,  with  the  00  array  being  the 
farthest  north.   The  Range  axis  is  191°  18'  1^.1"  T. 
2.   Computer  System 

The  computer  system,  a  Scientific  Data  Systems  (SDS) 
model  920,  is  installed  on  Zelatched  Point  and  consists  of 
three  main  subsystems:   the  signal  processing  subsystem, 
the  data  collection  subsystem,  and  the  computer  subsystem. 
The  signal  processing  subsystem,  the  link  between  the  hydro- 
phones and  the  data  collection  subsystem,  processes  the 
multiplexed  signal  received  from  the  arrays,  discriminates 
against  unwanted  signals,  and  determines  the  in-water  transit 
time  to  each  hydrophone.   The  data  collection  subsystem 
establishes  the  master  clock  timing  by  which  all  timing  is 
computed.   It  also  prepares  all  the  array  signals  for  the 
computer  subsystem,  which  calculates  the  tracked  object's 
position,  prints  the  tracking  data,  and  records  the  data  on 
magnetic  tape.   For  real-time  examination,  the  data  are 
plotted  on  X-Y  plotters. 

B.   EXPERIMENTAL  PROCEDURE 

A  92.5  dB  acoustic  source  with  approximately  ten  pounds 
of  ballast  was  suspended  from  the  stern  of  the  sound  boat 
to  25,  50,  75,  100,  150,  and  200-foot  depths  at  six 
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horizontal  ranges  (stations)  from  900  feet  to  4700  feet 
from  the  02  hydrophone  (Fig.  2).  The  first  station  was 
located  approximately  900  feet  north  of  hydrophone  array  02, 
and  450  feet  east  of  the  Range  axis.   Successive  stations 
were  located  parallel  to  the  axis  of  the  Range;  the  second 
being  directly  over  the  array,  the  third  600  feet  past  it 
(diametrically  opposite  the  first  station),  and  then  at  about 
1500-foot  intervals  to  about  4500  feet,  beyond  which  signal- 
to-noise  problems  made  further  examination  impossible.   The 
source  emitted  one  1.3  millesecond  75-KHz  pulse  per  second 
for  sixty  to  ninety  seconds  at  each  depth  at  each  station. 
The  position  of  the  boat  on  the  Range  was  found  in  the 
horizontal  from  the  Range's  Autotape  system,  a  portable  and 
commercial  microwave  system  used  exclusively  for  ship 
tracking  that  yields  accuracies  of  ±3  feet.   The  position 
of  the  acoustic  source  with  respect  to  the  Autotape  antenna 
was  determined  from  the  length  of  cable  let  out,  the 
horizontal  distance  from  the  antenna  to  the  stern  of  the  boat, 
and  the  boat's  heading.   Tidal  data  were  recorded  and  used 
for  correcting  depth  coordinates  to  a  mean  low-low  water 
(0-foot  tide)  reference  plane.   The  array  is  known  to  be 
585  feet  below  MLLW,  so  that  the  depth  of  the  source  below 
the  surface  and  the  tide  measurement  allows  calculation  of 
the  vertical  distance  between  the  source  and  array.   All 
measurements  were  made  with  respect  to  the  array's  coordinate 
system,  although  the  array's  orientation  and  location  on  the 
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Range  were  known,  making  it  possible  to  transform  this  data 
to  a  true  position  on  the  Range.   The  X  axis  of  the  02  array 
was  25°  24'  36. 0"  T. 

A  NAVTORPSTA  digital  velocimeter  that  has  an  accuracy  to 
±1.0  ft/sec  and  resolution  of  ±0.7  ft/sec  measured  the  sound- 
velocity  profile  at  each  station.   NAVTORPSTA  instrumentation 
at  Dabob  Bay  measured  and  recorded  the  acoustic  travel  times 
from  the  source  to  each  of  the  four  hydrophones  in  the  array 
to  six  decimal  places  for  every  source  depth  at  each  station. 
The  three-dimensional  position  of  the  source  with  respect  to 
the  array  was  calculated  from  these  travel  times  by  the 
NAVTORPSTA  NUTRACK  III  computer  program. 

The  data  were  gathered  on  21  March,  1974  during  a  neap 
tidal  period  and  while  the  sea  surface  was  smooth  and  glassy, 
although  there  had  been  ripples  earlier  in  the  morning.   The 
day  was  bright  and  sunny  with  scattered,  thin,  high  overcast, 
the  temperature  was  55°F,  and  no  noticeable  wind  was  present. 
The  shiny  stainless  steel  transducer  disappeared  from  view 
when  about  eight  feet  deep,  indicative  of  thick  biological 
concentrations  (gelbstaf fen)  in  the  water.   Rudimentary  drift 
measurements  of  boat  movement  relative  to  surface  water  were 
taken  by  throwing  paper  into  the  water  and  watching  its 
relative  motion  with  respect  to  the  boat,  which  indicated 
a  value  of  approximately  0.3  ft/sec. 
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III.   THEORETICAL  DEVELOPMENT 

The  computer  program  analysis  includes  corrections  for 
the  fact  that  acoustic  rays  are  not  straight  lines  in  a 
non-homogeneous  medium.   Two  simple  ray-tracing  techniques 
exist:   the  water  mass  divided  into  layers  of  isovelocity 
water  or  into  layers  with  a  constant  velocity  gradient. 
The  isovelocity  method  is  employed  in  the  NAVTORPSTA 
NUTRACK  III  program;  while  employment  of  the  isogradient 
technique  is  one  of  the  purposes  of  this  present  experiment. 

A.   NUTRACK  III 

1.   Raw  Array  Coordinates 

The  geometry  involved  in  the  calculations  of  the 
raw  array  coordinates  is  shown  in  Fig.  3.   As  described 
earlier,  the  four  hydrophones  are  on  four  adjacent  vertices 
separated  by  a  distance  d,  along  the  edge  of  the  cube.   The 
origin,  0,  of  the  array  coordinates  is  at  the  center  of  the 
cube  with  the  orthogonal  coordinates  parallel  to  its  edge. 
Let  the  transit  times  of  the  acoustic  pulse  from  the 
tracked  source  to  the  c  hydrophone  and  the  Y  hydrophone 
be  T   and  T  ,  respectively,  then 


(cTc)2  =  (RAWY  +  |)2  +  (RAWX  +  f)2  +  (RAWZ  +  |)2         (l) 

and 

(cT  )2  =  (RAWY  -  |)2  +  (RAWX  +  |)2  +  (RAWZ  +  |)2  .       (2) 
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2       2 
Taking  the  difference  between  T  and  T  ,  expanding, 

c      y 

and  combining  terms  yields 

2 
RAWY  =  §£  (Tc+Ty)(Tc-Ty)   .  (3) 

In  similar  fashion,  analogous  expressions  may  be 
obtained  for  the  X  and  Z  coordinates : 


c 


2 


and 


RAWX  =  —  (T  +T  )(T  -T  )  (4) 

2d   c    x   c    x  v  ' 


RAWZ  =  h   (Tc  +  Tz)(Tc-Tz)   >  (5) 


where 

c  =  the  average  sound  velcoity,  about  4860  ft/sec; 
d  =  the  distance  between  hydrophones,  30  ft; 
T  ,T  ,T  ,T  =  travel  times  in  seconds. 

The  travel  times  are  determined  from  a  clock  whose 

period  is  0.4x10"   sec.  (2.5  mHz),  which  is  sufficiently 

high  to  measure  precisely  the  arrival  time  of  a  75-kHz 

signal  (about  8  counts  per  quarter  period).   A  correction 

factor  must  be  applied  to  correct  for  the  array  receiver 

integration  (pulse  width  descrimination) .   Let 

I   =  the  integration  time,  800  us 

T'  =  the  measured  travel  time  to  the  R  hydrophone; 

T   =  actual  travel  time, 
c  ' 

so  that 


VI 


T»  =  T   +1   , 
c    c      ' 


Ty  "  Ty  +  T     - 


and 


Tc  +  Ty  =  (Tc  +  Ty  "  2I)   • 


T  —  T   =  T '  —  T ' 
c   y    c    y 


The  relations  for  the  raw  coordinates  are  now 

„2 
RAWY  =  Ig-  (T£  +  T£-K)(TJ~T£)   ,  (6) 


where  K  =   21   =   1.6x10    J   sec   and  T'    and  T'    are   the  measured 


"y 


t  ime  s . 
Similarly 


RAWX  =  %-?   (T1  +T»  -K)(T'  -T1)  (7) 

2d   c   y      ex 

and 

n2 

RAWZ  =  — ■  (T1  +T'  -K)(T'  -T')   .  (8) 

2d   c   z      c   z 


2 .   Tilt  Correction  to  Determine  Corrected 
Array  Coordinates 

The  raw  coordinates  previously  computed  are 

referenced  to  a  non-horizontal  reference  plane.   The  array 

is  instrumented  to  sense  the  angles  the  array  axes  make 

with  the  horizontal,  called  XTILT,  or  rotation  about  the  Y 

axis,  and  YTILT,  or  rotation  about  the  X  axis,  respectively 
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These  are  given  from 

XTILT  =  (TILTX  -  505)  x  .00833  degrees 
and 

YTILT  =  (TILTY  -  505)  x  .00833  degrees 

where  TILTX  and  TILTY  are  sensor  outputs  and  505  is  an 
arbitrarily  determined  sensor  reading  for  a  perfectly 
horizontal  orientation. 

The  correction  of  RAW  coordinates  for  array  tilt 
is  merely  a  matter  of  geometry.   Figure  4  shows  that 

CORX  =  [RAWX] [COS (XTILT)]  =  [RAWZ ] [SIN (XTILT) ] .  (9) 

Since  XTILT  is  small,  it  may  be  assumed  that  COS (XTILT)  =  1. 
The  relationship  then  becomes 

CORX  =  RAWX  -  [RAWZ] [SIN (XTILT)]   .  (10) 

Similarly 

CORY  =  RAWY  -  [RAWZ] [SIN (YTILT)]  (11) 

and 

CORZ  =  RAWZ  +  [RAWX][SIN(XTILT)]  +  [RAWY] [SIN (YTILT) ] 

(12) 

These  relationships  are  valid  for  both  positive  and  negative 

tilt  as  the  SIN  term  also  changes  sign  to  maintain  the 

relationship. 
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3.   Ray  Path  Refraction 

In  a  homogeneous  body  or  water,  the  minimum  time 
path  between  an  output  transducer  and  receiving  hydrophone 
is  a  straight  line.   When  the  speed  of  propagation  is 
spatially  varying  due  to  temperature  and  salinity  variations, 
refraction  occurs  and  the  minimum  time  path  is  no  longer  a 
straight  line.   Corrections  must  then  be  made  to  determine 
the  true  position  of  the  source. 

Prom  the  classical  laws  of  Fermat  and  Snell,  a 
minimum-time  path  is  computed,  working  back  in  time  from  the 
array  center  in  the  direction  of  the  "apparent"  position  and 
tracing  the  array  through  each  velocity  layer.   The  process 
is  continued  upward  in  space  and  backward  in  time  until  the 
measured  travel  time  is  consumed.   To  determine  ray-path 
refraction,  the  sound  speed  is  taken  at  one-foot  increments 
and  averaged  over  ten  feet.   Since  the  array  depth  is 
585  feet,  there  are  58  ten- foot  layers  and  one  five-foot 
layer.   The  vertical  angle  of  entry  of  the  ray  at  the 
array,  computed  from  the  tilt-corrected  array  coordinates 
(Pig.  5),  is 


A  =  SIN"1  [  ,-525*! ,  ]1/2  .      (13) 

0  CORX  +  CORY^  +  CORZ^ 


The  entry  angle  A   computed  from  the  CORR  array 
coordinates  is  assumed  to  be  equal  to  the  entry  angle  for 
the  minimum  time  path.   This  assumption  is  good  if  the 
source  is  far  enough  so  that  the  wave  fronts  are  plane  over 
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the  dimensions  of  the  array,  and  the  speed  of  sound  is 
constant  over  the  dimensions  of  the  array. 

The  travel  time,  T,  to  the  origin  is  computed  from 


Y"'fT  >  (l4) 

c    c 


where 


and 


S=  [CORX2  + CORY2  + CORZ2]1/2  (15) 


Rc  =  [(CORX +  15)2  +  (CORY  +  15)2  +  (CORZ +  15)2]1/2  , 

(16) 

where  R  is  the  slant  range  to  the  C  array  hydrophone, 
and  S  is  the  slant  range  to  the  array  center. 

It  may  be  seen  that  the  assumption  of  an  average 
velocity,  c,  in  computing  RAWX,  RAWY,  RAWZ  is  not  sensitive 
as  the  ratio  of  S  to  R  is  what  is  really  important  and  any 
small  error  will  tend  to  cancel  out. 

The  horizontal  distance  is 


H  =  [CORX2  +  CORY2]1/2   .  (17) 


The  first  layer  processed  Is  between  the  deepest 
depth  in  the  velocity  profile  that  is  shallower  than  the 
array,  and  the  array  itself; 

AZ1  =  Ad—  Z1         ,  (18) 

where  A,  is  the  array  depth. 
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The  incremental  slant  range  is 


AZ.. 

ASi  ■  side     '  (19) 


The  Incremental  travel  time  in  this  layer,  used  to  increment 
the  pulse  travel  time  T,  is 


AS 
AT  =  — -i   ,  (20) 

cl 


where  c-.  is  the  average  velocity  in  the  layer.   This 
incremental  travel  time  is  used  to  decrement  the  pulse 
travel  time  T.   The  horizontal  distance  for  this  time  is 


AH1  =  AS1   COS  AQ   .  (21) 


The  values  AH  and  AZ  are  accumulated  as  the  time  clock  is 
decremented  layer  by  layer. 

Using  Snell's  law,  the  change  in  angle  due  to 
refraction  upon  entry  into  a  new  velocity  layer  is 


COS  A      c 
___^ o_  _   1 

C0S(A^  +  AA)  '  c0 
o         2 


or 


c„ 
COS  A1    =  COS  A   —  (22) 

1  cl 


22 


where 


Al  =  Ao  +  AA 


and 


SIN  A1   =    (1  -  COS2  A1)1/2  (23) 


The  thickness  of  the  layer  is 


AZ2  =  Z2  -  Z]_    .  (24) 


This  process  is  repeated  and  the  original  travel 
time  T  is  decremented  with  the  incremental  layer  travel 
times  until  T  is  consumed.   The  ray-bent  X,  Y,  and  Z 
(RBNX,  RBNY,  RBNZ )  are 


RBNX  =  (  ^  )   Z  AH   ,  (25) 


RBNY  =  (  £^  )   Z  AH    ,  (26) 

and 

RBNZ  =  £  AZ  (27) 

This  raypath  migration  process  assumes  a  laterally 
homogeneous  medium  within  the  neighborhood  of  each  array. 
Independent  velocity  profiles  are  not  taken  at  each  array. 
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This  will  introduce  small  discontinuities  in  track  as  the 
tracked  vehicle  is  passed  from  one  array  to  another. 
4 .   Error  Estimation 

Returning  to  Eq.  A,  and  assuming  equal  errors, 
AT  =  AT   =  At   in  time  measurement:  the  maximum  expected 

CX  }  r 

error  in  RBNX  (AX)  is  found  from 


2 

RAWX  +  AX  =  %-?   [(T  +At)2-(T  -At)2] 
2d     c  x       J 


=  §5  ^T2-T2]  +  ^  [TcAt+TxAt]  (28) 


to  be 


AX^  [Tc+Tx]At   .  (29) 


At  large  ranges,  T   «  T  ,  so 

C       X 


2c 


AX  =  — -  T  At   , 
d   c 


but  c-T   =  R  ,   so 
c      5 


R      vd 
Similarly, 


AX  =  2  (#)  At   .  (30) 


and 


£±-  =  2  (§)  At  (3D 


^jr  *  2  (|)  At   .  (32) 


2H 


Assuming  the  maximum  error  in  time  is  one-fourth  of  the 
period  of  the  75-kHz  signal,   At  a  3.3x10   sec,   so 


AX    AY    AZ    n   nn-3 
IT  "  IT  m   IT   *    1  X  10 


B.   ISOGRADIENT  PROGRAM  "STUTRACK  I" 

An  isogradient  computer  program,  STUTRACK  I,  was 
developed  for  comparison  with  NUTRACK  III. 

1.   Linear  Gradient  Theory 

Figure  6  shows  a  highly  exaggerated  speed  of  sound 
profile  which  has  been  approximated  by  two  layers  each 
possessing  a  constant  gradient.   In  each  layer,  the  speed 
of  sound  may  be  given  as 


c(Z)  =  cQ  +  GAZ  (33) 


where  G  is  the  gradient;  G  is  negative  in  the  upper  layer 
and  positive  in  the  lower  layer.   Applying  Snell's  Law  for  A, 


C0S  A  =  CONSTANT  ,  (3*1) 


one  may  see  that  the  raypath  curves  downward  in  the  upper 
layer  and  upward  in  the  lower  layer.   Let  the  radius  of 
curvature  for  the  arc  between  the  angles  A,  and  A„  in  the 
upper  layer  be  R.   Then 
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C1"C2 


AZ   =   R(COS   A,  -COS    AP)    =    R  COS  A-     (-^ -)     .  (35) 

However,  from  Eq.  33, 


c1  -  c2   =   -GAZ 


and  with  Eq.    3H, 


R  -  t  coir    •  (36) 


Thus  for  each  raypath  the  value  of  R  is  a  constant 
in  a  linear  gradient  layer,  and  the  path  is  an  arc  of  a 
circle.   The  center  of  the  circle  lies  at  the  depth  for 
which  A  would  become  90°  if  the  sound  profile  for  the  upper 
layer  were  to  descend  to  greater  depths. 

Assuming  that  the  initial  angle  of  depression  of  the 
raypath  is  A  ,  from  the  geometry  of  Fig.  6  and  with  Eq.  36, 
the  change  in  depth  and  horizontal  range  become 


c 
AZ  =  -  Q  CQg  A   (COS  AQ  -  COS  A)  (37) 

and 

Ar  =  -  0  Cog  AQ  (SIN  A  -  SIN  V   >  (38) 


where  Z  is  the  increase  in  depth  and  r  is  the  increase  in 
range  from  the  reference  position,  and  A  is  the  angle  of 
depression  of  the  raypath  at  this  range  and  depth. 
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When  AZ  is  known,  the  new  angle  A  may  be  found  by 


A  =  COS-1  [COS  A  (^  +  1)]   .  (39) 

o 


The  time-of-flight  along  the  raypath  may  be 
simplified  to 

A 

Ao 
which  upon  integration  yields 


..         1    ,    r,      COS   A      wl   +    SIN   Aon  ,IinN 

At   =  g  ln[(1  +   SIN  A)(      cos   A      ")]      .  (41) 


2.   STUTRACK  I  Construction 

STUTRACK  I  (App.  A)  is  designed  to  allow  the  deter- 
mination of  a  minimum  amount  of  environmental  information 
necessary  for  agreement  with  the  NUTRACK  III  results.   The 
program  determines  the  tilt-corrected  coordinates  and  the 
raypath  entry  angle  into  the  array  In  the  same  manner  as 
NUTRACK  III,  but  uses  Eq.  39  to  calculate  the  exit  angle  A 
from  the  layer,  Eq.  HO   to  calculate  the  horizontal  range 
traversed,  and  the  time-of-flight  from  Eq.  4l.   The  thickness 
of  the  layer,  AZ,  is  determined  on  a  point-to-point  basis 
from  the  depths  from  which  the  velocity  profile  was  established. 
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The  values  Ar  and  AZ  are  accumulated  and  the  original  travel 
time  T  Is  decremented  as  In  NUTRACK  III  until 

TG  >  T  +  e     , 


where 


and 


TG  =  Z  At 


e  =  .0005  sec. 

This  value  of  e  ensures  an  accuracy  of  Z  AZ  and  Z  Ar 
within  ±3  inches.   If  TG  is  within  e  of  T,  the  process 
stops  and  the  X3  Y,  and  Z  coordinates  for  this  isogradient 
technique  (XG,YG,ZG)  are  calculated  as  in  NUTRACK  III, 


XG  =  ^^  Z  Ar   ,  (42) 


YG  =  ^pl  Z  Ar   ,  (43) 


and 

ZG  =  ^P-   Z  Ar   .  (44) 


When 


TG  >  T  +  e    , 

the  program  returns  to  the  previous  layer  and  a  bracket- 
and-halving  process  begins  in  AZ  until  the  proper  AZ  is 
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found  to  satisfy  the  e  criterion  in  time.   XG,  YG,  and  ZG 
are  then  found  by  Eqs.  42,- #3,  \md'"«y... 

STUTRACK  I  is  also  designed  to  permit  the  utilization 
of  the  velocity-profile  data  in  varying  incremental  step 
sizes  to  allow  a  determination  of  the  effect  of  a  varying 
amount  of  environmental  data.   For  example,  if  velocimeter 
data  is  gathered  in  two-foot  steps,  an  incremental  step  size 
(INC)  of  1  would  automatically  set  AZ  =  2  ft,  and  an  INC  of 
10  would  set  AZ  =  20  ft.   Consideration  of  the  source  depth 
must  be  made  prior  to  the  invocation  of  STUTRACK  I  in  its 
present  form,  as  the  setting  of  too  large  an  INC  can  create 
I  AZ  greater  than  the  array  depth  of  585  feet,  resulting  in 
an  error  message. 

STUTRACK  I  is  also  designed  to  analyze  the  data  with 
the  same  isovelocity  technique  as  NUTRACK  III,  but  in 
varying  incremental  step  sizes  for  comparison  with  the 
isogradient  method.   These  output  coordinates  are  XI,  YI, 
and  ZI. 

The  inputs  into  STUTRACK  I  are:   velocity  profile, 

increment  size,  tilt  correction,  time  of  day  in  hours, 

minutes   and  seconds,  T  ,  T  ,  T  ,  T  ,  RBNX,  RBNY,  and  RBNZ. 

x   y    z    c 

The  output  Is  six  histograms  of  the  differences  between  the 
NUTRACK  III,  ray  bent  coordinates,  and  their  respective 
values  as  determined  by  STUTRACK  I  using  both  the  isogradient 
and  isovelocity  methods.   These  differences  and  the  STUTRACK  I 
values  are  also  printed. 
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IV.   EXPERIMENTAL  RESULTS 

A.  VELOCITY  PROFILE 

Figures  7  and  8  Illustrate  that  the  velocity  profiles 
for  station  1  (taken  at  1115  hrs)  and  station  6  (taken  at 
1315  hrs)  are  equivalent  below  85  feet.   (See  Fig.  2  for 
station  positions.)   The  data  for  station  6,  represented  by 
the  "+"'33  shows  the  "afternoon  effect"  of  solar  heating. 
Both  profiles  increase  in  velocity  with  depth  at  a  greater 
rate  than  would  be  caused  by  pressure  alone,  implying  a 
similarly  increasing  temperature  profile.   The  colder 
surface  water  was  fresher,  however,  due  to  the  heavy  spring 
rainfall  and  runoff  from  the  snow-capped  mountain  peaks. 

B.  AUTOTAPE  DATA 

The  Range's  microwave  Autotape  system  was  used  to  track 
the  sound  boat.   Its  precision  was  much  less  than  expected 
so  that  comparison  with  Range  data  was  not:  attempted  at  this 
time  and  will  have  to  wait  until  the  Autoftape  data  is 
subjected  to  extensive  analysis. 

C.  ANALYSIS  OF  NUTRACK  III  RESULTS 
1.   One-Way  Transmission  Loss 

The  one-way  transmission  loss  (TL,)  is  determined  by 
an  assumption  of  spherical  spreading  and  "the  relationship 
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TL(dB)  =  20  LOG  R  +   aR     ,  (45) 

where 

a  =  attenuation  constant  in  sea  water 
for  75-kHz  sound, 

-    0.03  dB/meter. 

Figure  9  shows  that  at  short  ranges  (station  1)  the  TL  is 

dominated  by  the  spherical  spreading  term.   Beyond  a  "critical 

range"  of  approximately  1000  feet  for  75-kHz  sound,  the 

losses  due  to  absorption  of  sound  in  sea  water  will  overpower 

the  spreading  and  the  transmission  loss  will  increase 

rapidly  at  greater  distances.   The  TL  from  station  6  to  the 

hydrophone  array  is  approximately  105  dB,  but  is  only  57  dB 

from  station  1. 

2.   RBNX,  RBNY,  RBNZ  and  Slant  Ranges 
For  Stations  1-25  and  6-"50 

Figures  10-13  show  the  80  data  for  station  1-25 

plotted  versus  time.   Figures  14-16  represent  the  same  data 

as  Figs.  10-12  but  with  a  linear  drift  correction.   Figures 

17-21  are  similar  curves  for  the  66   data  for  station  6-50, 

except  that  a  drift  correction  was  applied  only  to  the  slant 

range  curve.   These  two  stations  were  chosen  for  preliminary 

investigation  since  station  1  was  closest  to  the  array  while 

station  6  was  farthest  away.   The  RBNX,  RBNY,  and  Slant  Range 

scales  for  each  station  are  nearly  equivalent,  but  RBNZ  data 

are  plotted  in  a  5/1  ratio,  since  station  6  was  five  times 

the  slant  range  of  station  1.   In  all  cases,  station  1  data 
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show  good  point-to-point  continuity  while  station  6  data, 
except  for  Slant  Range,  are  very  erratic.   Microstructure 
was  examined  as  the  cause  of  this  variation  using  the 
procedure  of  Medwin  [2],  but  was  determined  to  be  capable 
of  producing  a  maximum  fluctuation  of  only  ±1  foot,  or  only 
a  small  fraction  of  the  actual  variation. 

3.   Transformation  to  Spherical  Coordinates 

A  transformation  from  the  orthogonal  coordinates 
RBNX,  RBNY,  and  RBNZ  to  a  spherical  system  was  achieved 
using  the  following  relations, 


R  =  (RBNX2  +  RBNY2  +  RBNZ2)1/2   ,  (46) 


$  =  SiN-1(^Z)   , 
n 


and 


u    iHN   ^RBNX; 


5 


where  $  is  the  angle  of  elevation  looking  from  the  array  to 
the  source,  and  0  is  an  arbitrarily  defined  azimuthal  angle. 
A  graph  of  <J>  versus  0  is  a  scatter  diagram  of  the  projection 
of  the  data  onto  the  surface  of  a  sphere  of  radius  R. 

Figures  22-29  are  station  1-25  plots  of  $  and  0  as 
calculated  from: 

a.  Ray  bent  coordinates  (Figs.  22-24), 

b.  Ray  bent  coordinates,  amplified,  and  with  linear 
drift  correction  (Figs.  25-27), 

c.  COR  coordinates,  amplified  (Figs.  28  and  29). 
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Figures  22-24  show  a  consistent,  nearly  linear  drift,  with 
good  point-to-point  continuity.   Figures  26  and  29  are 
identical,  as  should  be  the  case  since  NUTRACK  III  assumed 
a  laterally  homogeneous  medium,  thereby  keeping  the  ratio 
Y/X  constant.   Figures  2  5  and  2  8  show  similar  but  not 
identical  variations,  which  was  also  expected  since  the  COR 
coordinates  assumed  a  vertically  homogeneous  medium  and  the 
raybent  coordinates  consider  the  inhomogeneities. 

Figures  30-32  are  station  6-50  plots  of  $  vs.  time, 
0  vs.  time,  and  $  vs .  0,  respectively,  as  calculated  directly 
from  the  raybent  coordinates,  and  plotted  to  the  same  scales 
as  Figs.  25-27.   The  variations  are  approximately  the  same 
in  magnitude  as  for  the  same  plots  at  station  1-25.   However 
the  scatter  diagram  for  station  1-25  transforms  into  a  five- 
foot  square  patch  on  the  surface  of  a  sphere  of  radius  900 
feet.   At  station  6-50  this  uncertainty  grows  with  the  ratio 
of  the  radii  to  a  patch  twenty-five  feet  square. 
k .   Error  Determination 

An  error  determination  was  made  from  the  data  from 
Figs.  12,  1H,  15,  and  17-19  by  subtracting  each  datum  from 
the  average  value  drift-corrected  coordinate  and  dividing 

by  R.   Variations  from  the  averages  were  found  to  be  less 

—  6  —  3 

than  5  x 10"   with  a  standard  deviation  (a)  less  than  10 

This  latter  is  comparable  with  the  theoretical  error  deter- 
mination of  section  III-A-4  on  the  basis  of  an  error  in 
determining  the  arrival  time  which  is  independent  of  range. 
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D.   COMPARISON  OF  NUTRACK  III  TO  STUTRACK  I 

Comparison  of  STUTRACK  I  data  in  both  the  isogradient 
and  isovelocity  modes  with  the  NUTRACK  III  raybent  coordinates 
was  achieved  by  the  following  relationships 

DXG  =  RBNX  -  XG  (47) 

DYG  =  RBNY  -  YG 
DZG  =  RBNZ  -  ZG 
and 

DXI  =  RBNX  -  XI  (48) 

DYI  =  RBNY  -  YI 
DZI  =  RBNZ  -  ZI   . 

Station  6-50  data  (Figs.  36-38)  are  consistent  with 
station  1-25  data  (Figs.  33-35),  but  much  more  dramatic 
since  it  was  taken  at  a  much  greater  range.   These  graphs 
are  plots  of  DX,  DY,  and  DZ  for  both  the  isogradient  and 
isovelocity  cases  as  a  function  of  incremental  step  size. 
They  show  that  the  STUTRACK  I  results  are  essentially 
equivalent  to  the  NUTRACK  III  data  when  INC  =  1.   However 
as  INC  is  increased,  the  isovelocity  technique  gives 
consistently  poorer  results  while  the  isogradient  method's 
results  remain  adequate  as  INC  increases.   In  all  cases  the 
STUTRACK  I  coordinates  followed  the  NUTRACK  III  raybent 
coordinates  point-for-point ,   but  were  biased  to  one  side. 
A  few  seconds  of  computer  time  was  also  saved  as  a  result 
of  the  isogradient  technique  at  a  greater  INC  than  for  the 
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comparable  calculations  via  the  isovelocity  method  since 
fewer  layers  had  to  be  examiTCfl  -to  "Bcnieve  the  same,  or 
better,  results. 

Figure  39  is  a  histogram  of  DXG,  INC  =  1 ,  at  station 
6-50.   Its  Gaussian  nature  is  characteristic  of  all  the 
histograms  of  DX,  DY,  and  DZ  at  every  INC  tested,  at  both 
stations. 
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V.   CONCLUSIONS 

A.  NUTRACK  III  DATA 

The  large  variations  in  RBNX,  RBNY,  and  RBNZ  compared  to 

the  smoothness  of  the  Slant  Range  may  be  explained  by  an 

examination  of  Eqs.  3,  4,  5,  and  15.   The  individual 

coordinates  are  dependent  upon  the  difference  of  T  and  T  , 

T„,  or  T  .   Since  this  difference  is  small,  a  small  error 
y     z  ' 

in  time  will  make  a  large  fractional  error  in  the  computed 
coordinate.   The  uncertainty  in  signal  discrimination, 
addressed  in  section  III-A-4  and  estimated  as  3-3x10"  sec, 
becomes  important,  causing  significant  random  errors  from 
pulse-to-pulse.   This  error  becomes  very  large  for  any 
position  for  which  the  wave  arrives  at  the  two  hydrophones 
at  nearly  the  same  time.   For  calculation  of  the  Slant  Range 
the  arrival  time  is  used  alone  so  the  fractional  error  is 
much  less. 

B.  COMPARISON  OP  STATION  1-25  TO  6-50 

Comparison  of  station  1-25  data  to  that  from  6-50  shows 
the  latter  to  have  greater  spatial  variability.   The  trans- 
formation to  spherical  coordinates  shows  this  difference  is 
not  present  in  the  angular  display.   This  is  what  would  be 
expected  if  the  time  errors  were  the  important  sources  of 
variability.   This  conclusion  is  upheld  by  an  experimental 
error  determination  at  both  stations  that  agrees  with  the 
theoretical  error  Induced  by  the  one-fourth  period 
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uncertainty  in  the  signal  arrival  times.   Out  to  the 
greatest  distances  studied  (over  twice  that  used  in  normal 
Range  operations),  the  precision  in  measuring  arrival  times 
is  unchanged  even  though  the  signal-to-noise  ratio  decreases 
by  nearly  50  dB.   Therefore,  it  is  concluded  that  for 
conditions  similar  to  those  prevailing  at  the  time  of  this 
experiment  any  attempt  to  improve  the  signal-to-noise  level 
by  increasing  the  source  level  will  be  nonproductive  and 
wasteful.   Improved  Range  accuracy  can  only  be  obtained  by 
improving  the  accuracy  of  the  time  measurements  (above  its 
present  value  of  about  one-quarter  period  of  the  75-kHz 
signal)  or  increasing  the  baseline  of  the  hydrophone  array. 

C.   EFFECT  OF  ENVIRONMENTAL  DATA 

The  graphs  of  DX,  DY,  and  DZ  for  both  stations  show 
conclusively  that  for  the  conditions  prevailing  at  the  time 
of  this  experiment  the  isogradient  technique  of  STUTRACK  I 
provides  better  results  with  less  environmental  data  than 
does  NUTRACK  III.   This  is  consistent  in  each  of  the  plots 
but. is  explicit  in  DZ  at  station  6,  where  the  isovelocity 
layer  method  results  in  a  linear  growth  in  DZ  as  the 
incremental  step  size  (INC)  is  increased,  but  little  change 
in  DZ  in  the  isogradient  technique.   Conclusive  in  this  is 
that  much  less  environmental  data  need  be  collected  if  a 
isogradient  program  were  used. 
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VI.   COMMENTS  ggfg  -j^vSTB^ENDATIONS 

This  thesis  treated  the  NUTRACK  III  data  as  a  comparator 
for  analysis  for  both  modes  of  STUTRACK  I,  but  no  conclusive 
evidence  has  been  gathered  to  determine  which  is  more 
accurate.   The  Range's  Autotape  system  was  to  have  determined 
this,  but  was  so  imprecise  that  a  decision  was  unreachable 
without  detailed  analysis  of  the  Autotape  results. 

In  all  Z  calculations  at  both  stations  the  source  depth 
plus  Z  was  less  than  the  array  depth  by  three-to-ten  feet. 
No  hypothesis  has  been  reached  for  the  cause  of  the 
disappearance  of  these  few  feet. 

Further  study  of  the  data  is  needed  to  determine  if  the 
results  from  the  intermediate  stations  are  consistent  with 
the  conclusions  reached  here.   Then,  these  conclusions  should 
be  tested  using  sound-velocity  profiles  representing  other 
Dabob  Bay  oceanographic  conditions.   Special  consideration 
should  be  given  to  the  apparent  inaccuracy  in  Z. 


•50 


FIGURE  1.   DABOB  BAY  RANGE 
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FIGURE    5.       RAY    PATH    REFRACTION 
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FIGURE  6.   ISOGRADIENT  THEORY 
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FIGURE  7.   VELOCITY  PROFILE 
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FIGURE  8.   VELOCITY  PROFILE 
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FIGURE    9.       TRANSMISSION   LOSS 
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FIGURE    39.       DXG   HISTOGRAM,    Station    6-50 
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APPENDIX  A 
STUTRACK  I  COMPUTER  PROGRAM 
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